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ABSTRACT 

We use a multi-dimensional Monte Carlo code to compute X-ray spectra for a variety of 
active galactic nucleus (AGN) disk-wind outflow geometries. We focus on the formation of 
blue-shifted absorption features in the Fe K band and show that line features similar to those 
which have been reported in observations are often produced for lines-of-sight through disk- 
wind geometries. We also discuss the formation of other spectral features in highly ionized 
outflows. In particular we show that, for sufficiently high wind densities, moderately strong 
Fe K emission lines can form and that electron scattering in the flow may cause these lines 
to develop extended red wings. We illustrate the potential relevance of such models to the 
interpretation of real X-ray data by comparison with observations of a well-known AGN, 
Mrk 766. 
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1 INTRODUCTION 

The study of active galactic nuclei (AGN) is an important topic in 
contemporary astrophysics. These objects are interesting in their 
own right, allowing us to study the extreme physics of accretion 
in the vicinity of a supermassive black hole. Furthermore, numer- 
ical simulations suggest that AGN likely have a critical role in the 
formation and evolution of galaxies, highlighting the need to under- 
stand how these objects accrete matter, grow and feedback energy 
to their surroundings (e.g. ,Croto n et al. 2006). 

Since AGN are bright X-ray sources, they have been pop- 
ular targets for almost all X-ray observatories. One of the most 
interesting results obtained thanks to the high sensitivity of the 
current generation of X-ray m issions (specifically X MM-Newton 



Ijansen et alj|200l]l. Chandra 



recently, Suzaku iMitsuda et al 



Weisskopf et al. 
l2007ll ) 

narrow absorption features in the 2-10 keV band of several br ight 
AGN (see, e.g.jPou nds et al. 2003 ; Reev es et a l. 2004 ; Risaliti et^ 
l2005l : lYoung etafi l2005l : [Turner et al.l |2007j; iBraito et alj l2007h . 
Particularly for features in the Fe K region of the spectrum, the ob- 
served energies and strengths of these lines suggest identification 
with very highly ionized species (e.g. Fe XXV and XXVI) in very 
fast (up to ~ 0.1c) outflows. Blueshifted absorption lines are well- 
known from observations of AGN in other wavebands and models 
explaining suc h phenomena in terms of winds have been devel- 
oped (see e.g. iMurrav et ai]|l995l lElvisI |2000|. IProga & KallmaiJ 
|2004) . However, the new X-ray data clearly suggest a compo- 



nent of very highly ionized fast outflowing plasma, the relation- 
ship of which to the less ionized material observable in ultravio- 
let (uv) or softer X-ray wavebands remains unclear. In particular, 
moderately ionized outflows identified in soft X-ray spectra gener- 
ally have somewhat lower ve locity ( < 1000 km s~^: ,Blustin et al] 
I2OO5I iMcKernan et alJbOOTh while uv spectr a can show evidence 
of either low or high velocity absorption (e.g. lElvisll2000l) . 

As observational evidence in support of the phenomenon has 
accumulated, several theoretical studies of the physi cal properties 
of highly ionized AGN outflows have been made. IPounds et al.l 
(2003) and King & Pounds (2003) discussed the blueshifted ab- 
sorption features in PG 12 11+143 in terms of a conical outflow 
subtending a large solid angle. They suggested that such a flow 
might be driven by continuum radiation pressure and that emis- 
sion from such a flow might be responsible for the big blue bump 
which dominates the bolom etric output of PG121 1+143. However, 
lEverett & Ballantynd j2004l) considered flows driven by continuum 
radiation pressure in greater detail and concluded that, while this 
mechanism could work in principle, the resulting outflows were 
un likely to produce sp ectral signatures as strong as those reported 
bv |Poundsetai] ( l2003h . 

[siii] ( l200^ undertook a two-dimensional (2D) Monte Carlo 
(MC) radiative transfer study of parameterised conical outflow 
models and concluded that, for suitable column densities, viewing 
down such flows could account for the observed blueshifted, nar- 
row absorption lines in PG121 1 + 143. However, that study was lim- 
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ited to consideration of only the simplest conical geometry and did 
not consider the effect of either rotation or off-axis lines-of-sight 
on the spectrum. 

Using a chained ve rsi on of the ID XSTA R code 
jKallman & Bautistal 1200 ih , ISchurch & Don3 ( l200l l2008h 
have computed spectra for columns of outflowing gas with a 
variety of density and velocity profiles. They demonstrated that 
outflows can imprint a wide variety of features on X-ray spectra, 
depending on the outflow conditions and that very high outflow 
velocities would be required if absorption in outflows is to explain 
soft excesses in A GN spectra. More recen tly, a similar approach 
to that adopted by ISchurch & Pond ( l2007h has been used for the 
calculation of transmission spectra i n dynamical models of AGN 
outflows bv lPorodnitsvn et al.l ( l2008h . 

In this paper, we extend the study of ISimI to mcorpo- 

rate a more realistic and more versatile description of a disk wind 
geometry including both rotation and off-axis lines-of-sight. We re- 
tain the use of the MC method owing to its versatility for multi- 
dimen sional radiative transfer . Our method is complementary to 
that of is church & Don3 ( l2007h since we account for geometric ef- 
fects directly (e.g. scattering of radiation between lines-of-sight in 
multi-dimensional outflow geometries) but make some simplifica- 
tions in the treatment of atomic processes. 

We focus on highly ionized winds and the interpretation of ob- 
servable features in the Fe K band. Although other spectral lines are 
expected to form in outflows, Fe Ka absorption is the most relevant 
for the interpretation of observed spectra since it is the clearest sig- 
nature of a highly ionized flow: although more abundant, the lighter 
elements are fully ionized more easily than Fe so that any spec- 
tral features they imprint are weaker. Furthermore, identification of 
features in the Fe K energy band is relatively secure since only K 
shell transitions of heavy ions are expected at energies ^ 5 keV. 
At lower energies, where Ka transitions of light or intermediate 
mass elements might arise, line identification is more ambiguous, 
particularly if large velocity shifts are considered. 

Although we regard Fe Ka absorption features as the best 
diagnostic for a highly ionized flow, we will also investigate the 
formation of emission components in the Fe K arising from line 
scattering or recombination in the outflow. Such features are of 
particular interest sin ce th ey m ay e xplain possible P- Cygni-like 
line profiles (e.g. IPone et al.ll2007l . iTurner et alj|2008h and may 
affect the interpretation of Fe emissio n features commonly asso- 
ciated with disk reflection (see e.g. iLaming & TitarchuM 1 20041 
iLaurent & Titarchukll 20071) . 

We begin in Section[2]by defining the class of outflow mod- 
els which we will consider. In Section[3]we describe our radiative 
transfer method and discuss the adopted atomic data. We present 
a sample calculation for one model in Section |4] and then extend 
our discussion to a a grid of outflow models in Section[5] The im- 
plications of our models for observations of Fe Ka absorption are 
described in Section|6]and for other spectral features in Section|7] 
In Section [8] we illustrate the value of our models by comparing 
them in detail with observations of a well-known AGN, Mrk 766. 
We draw conclusions and discuss further work in Section|9] 



2 MODEL 

Our radiative transfer calculations were performed using a simply- 
parameterised model for an outflow launched from an accretion 
disk around a supermassive black hole. In this section, we describe 
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Figure 1. The geometrical construction used to define the wind (only the 
positive xz-plane is shown - the wind is symmetric under both reflection 
in the a;j/-plane and rotation about the z-axis). The region occupied by the 
wind is shaded. The three parameters which determine the geometry {d, 
^min, ''max) are indicated. The star symbol represents the point in the wind 
which is specified by the wind coordinates R and I. 

the properties of the model and the parameters which must be spec- 
ified to define a particular realisation of the model. 

2.1 Geometry 

We adopt a stand ard disk wind geome try, namely the "displaced 
dipole" model of iKnigge et al.l ( Il995h . This geometry has been 
adopted in radiative transfer studies of accretion dis k winds for a 
variety of systems inc l uding cataclysmic variables jKnigge et al.l 
1995; Long & Kniggd l2002h and massive young stellar objects 
( Sim et al.,,2005i) . The geometry is illustrated in Fig. [T] and is de- 
fined by the following three parameters (each of which is marked 
in the figure): 

(i) d, the distance of the focus point below the origin 

(ii) J'min, the distance from the origin to the inner edge of the 
wind in the xy-plane 

(iii) '"max, the distance from the origin to the outer edge of the 
wind in the xy-plane 

The accretion disk is assumed to lie in the xy-plane. The wind 
is symmetric under rotation about the z-axis and under reflection in 
the a;i/-plane. 

2.2 Velocity law in the vtind 

The velocity is specified a t every point in the wi nd following the 
param eterisation of Knig ge et al.l l ll995h (see also lLong & Kniggd 
I2OO2I) . Although this velocity law is not based on self-consis tent 
hydrodynamical outflow models (cf . iDorodnitsvn et al.l 1200 Sh . it 
provides a simple and reasonably flexible description of possible 
steady-state flows which we use for our exploratory radiative trans- 
fer simulations. 

2.2.1 Rotation 

The wind rotates about the z-axis. The rotational velocity is ob- 
tained by assuming that parcels of matter conserve specific angular 
momentum about the z-axis as they flow outwards. The angular 
momentum at the base of a streamline is set at the Keplerian value 



3 



for the radius at which the streamline crosses the xy-plane. Thus 
the rotational velocity is determined only by the choice of wind 
geometry (see above) and the mass of the central object, Mhh- 



2.2.2 Outflow 

The outflow velocity points directly away from the focus point of 
the wind. Its magnitude is given by 



Vl — Vo + [Vo 
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where / is the distance along the outflow streamline (see Fig. [T](, 
Rv is an acceleration length parameter and the exponent, /3, sets 
the rate of acceleration. The terminal velocity, v^o , is specified as a 
multiple /„ of the escape speed from the base of the streamline 
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where R is the radius at the base of the streamline (see Fig.lTJ. 



2.5 Atomic level populations 

In order to keep the computational requirements of our multi- 
dimensional MC calculations tractable, we have not attempted a 
complete non-LTE calculation, but have instead employed approx- 
imate formulae for ionization and excitation. Our approach is sim- 
ilar to the standard modified nebular approximation, which has 
been used in a variety of ot her astroph ysical calculations (e.g. stel- 
lar w inds lAbbott& Lijcvl ll985. Luc y & Abbot^l 19931 Fvirik et al.l 
Il999ll . supemovae fMazzali & Lucv 1993] and disk winds from 
cataclysmic variables [Long & Knigge 2002]). However, as dis- 
cussed in Section [J!4l we have modified standard approximations to 
reflect the fact that X-ray spectra of AGN are better approximated 
by a power law than a blackbody. 



2.5.1 Ionization 

For simplicity, we assume ionization equilibrium in which the dom- 
inant ionization process for all the ions we consider is photoioniza- 
tion from the ground state. Furthermore, we neglect recombinations 
from excited states such that the condition of ionization equilibrium 
can be expressed: 



2.3 Mass density 

The wind is assumed to be smooth and in a steady state flow; exten- 
sion of this study to clumpy flows will be the topic of a subsequent 
investigation. The total mass-loss rate (M) of the wind is treated 
as a parameter and the mass loading of specific outflow streamlines 
is determined by parameterising the mass-loss rate per unit surface 
area as a function of R via 



dm k 
— — oc R 
dA 



(3) 



subject to the constraint that 



(5) 



where riifi is the atomic level population of an ionization state i in 
its ground state (denoted 0), Ue is the free electron density, 7^,0 is 
the photoionization rate coefficient from the ground state of ion i 
and cti+ifi^i is the recombination rate coefficient from the ground 
state of ionization state i + 1 to level I of ionization state i. The 
summation runs over all states of the lower ion. 
This can be rewritten as 



(6) 
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dm 
U 



RdR = M 



(4) 



where $o.o,i(Te) is the population ratio ni+l,o?^e/7^i,o evaluated 
for LTE at the electron temperature (Te) and 



The mass density at a point in the wind is then obtained by com- 
bining the specific mass-loss rate on local outflow streamlines with 
the outflow velocity, assuming a stationary flow. 



(7) 



is the fraction of recombinations that go directly to the ground state. 



2.4 Electron temperature 

The electron temperature Te of the outflowing material must be 
specified since it has a role to play in determining the ioniza- 
tion/excitation state of the plasma (see below). Ultimately, it should 
be calculated by consideration of all relevant heating and cool- 
ing processes. But the treatment of all such processes - heating 
and cooling via lines, bound-free continua, free-free processes and 
(inverse-)Compton scattering of photons at all energies - would re- 
quire more detailed atomic physics, simulations of a much wider 
range of photon energies and the self-consistent consideration of 
heating by radiation from other sources (e.g. radiation emitted or 
reflected by the accretion disk or surrounding stru ctures). Suc h is- 
sues go beyond the scope of this study and so, as in lSimI j2005l) . Te 
is assumed to be uniform and is treated as an input parameter. 



Si{Te,J) = ^''lj^ = / duJuf ^ du/ I ayBv{Te)u ^ du 

(8) 

is the ratio of the true photoionization rate determined from the 
mean intensity (J^) and that which would be obtained in LTE at 
Te. ai, is the ground state photoionization cross-section for fre- 
quency I/, Bv is the Planck function and the integrations run from 
the threshold frequency, vo, to infinity. 

Equation [6] is used for all the ionization calculations in this 
paper. Since 'l>o,o,i and C, depend only on the adopted electron tem- 
perature these can be evaluated prior to the MC radiative transfer 
simulation. Si , however, depends on the radiation field Ju and so it 
is computed from the radiation field parameterisation and iterated 
(see Section [34t . 
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2.5.2 Excitation 

Since the majority of the atomic processes with which we are con- 
cerned are associated with ground state absorption, the treatment 
of excited level populations is relatively unimportant and there- 
fore very approximate. For all excited states, we adopt a two-level 
radiation-dominated Sobolev approximation 

I^^ELf^ldi + iY' (9) 
no go \ JuC^ J 

where hui is the excitation energy of state I above the ground state, 
gi is the statistical weight of level I, [3 is the angle-averaged Sobolev 
escape probability for the transition from level I to the ground state 
and Jv is the mean intensity in the transition. 



packets can interact with the outflowing plasma, thereby changing 
both their direction of propagation and their observer-frame fre- 
quency. In all interactions, energy conservation is strictly enforced 
such that, in a local co-moving frame, there is no net gain or loss of 
radiative energy. The particular types of packet interaction which 
may occur during the simulation are summarised below. 



3.3.1 Compton scattering 

By far the most common interaction is Compton scattering b y free 
electro ns. This process is treated using the method outlined bv lLucvl 
( l2005h . In all our simulations, we assume the electron temperature 
is small enough that inverse Compton scattering may be neglected 
in the co-moving frame. 



3 RADIATION TRANSPORT SIMULATIONS 



The propagation of radiation through the model is followed via a 
MC simulation in which the quanta are indivisible packets of ra- 
diative energy. The principle s of the method have been developed 



and described elsewhere (e.g. 


Abbott & Lucvl 19851 Lucv & AbbottI 


ll993llMazzaU & Lucvll 19931 


Lucvll2002llLucvl 2003b so onlv those 



points of specific relevance to the procedure used in this study are 
given below. 



3.1 Computational grid 

For the radiative transfer calculations, the wind properties are first 
discretised onto a grid in the natural wind parameters [R, I). Typ- 
ically a 100 X 100 grid is used. For convenience when propa- 
gating the packets, this grid of wind properties is then mapped 
to an underlying three-dimensional (3D) Cartesian grid, typically 
100x100x100. Each Cartesian cell which lies inside the wind is 
assigned the density, ionization/excitation state and radiation field 
properties of the closest point in the grid of wind properties while 
those Cartesian grid cells which lie outside the wind are assumed 
to be empty. 

3.2 Initialisation of packets 

Since our primary objective is to simulate spectra in the keV energy 
range, we initially create packets between 0. 1 and 40 keV. A power- 
law distribution of packet energies is assumed such that the input 
spectrum n{E) varies with photon energy (E) according to 

n(S)[photons s~^ keV"^] oc . (10) 

This is an adequate first-order description of the observed X-ray 
spectra of AGN in the 2-10 keV range with which we are primarily 
concerned. The power-law index (F) is as an input parameter for 
our model. 

It is assumed that the primary X-ray source is compact and 
roughly appears as a point source as seen from the outflowing gas. 
Thus, the packets are launched from the coordinate origin. Their 
initial directions are chosen isotropically. 

3.3 Propagation of packets 

Once launched, the packets propagate until they reach the outer 
boundary of the computational domain. As they propagate, the 



3.3.2 Photoabsorption 

Bound-free continua can absorb any photons with frequency above 
their edges - the cross-sections are obtained from the atomic data 
sources described in Section [33] During the MC simulations, only 
bound-free absorption by ionic ground states is included. 



3.3.3 Bound-bound absorption 

Atomic lines can absorb photons that come into resonance with 
them. Their optical depths are computed in the Sobolev approxi- 
mation using the velocity gradient and the appropriate level popu- 
lations obtained from the ionization/excitation formulae. 



3.3.4 Re-emission 

Following photoabsorpti on or bound-bound absorption, the macro 
atom scheme devised bv lLucvl ( I2002ll2003t) is used to determine 
the subsequent re-emission frequency of the packet. This approach 
naturally incorporates both resonance scattering and line emission 
following recombination, the two processes which are of greatest 
interest for outflow features in AGN X-ray spectra. The particular 
macro atom scheme adopted allows for bound-bound radiative and 
electron collisional processes in both internal state changes and de- 
activation. Bound-free recombination processes are also included 
but, in accordance with the assumption made for Equation|6] pho- 
toionization from excited states is neglected. 

Packets emitted by macro atoms via bound-free continua are 
radiated isotropically. For bound-bound transitions, the direction 
of re-emission is determined by constructing a 2D grid of Sobolev 
escape probabilities as a function of polar and azimuthal angle and 
sampling this distribution to choose a new packet trajectory. 



3.3.5 k-packets 

When the outcome of an interaction is a conversion of a radiative 
packet to a packets o f thermal energy (fc-packets in the nomen- 
clature of lLucvll2003h . we assume that the subsequent elimination 
of the fc-packets predominantly leads to emission at lower photon 
energies such that these packets are lost to the hard X-ray energy 
band. 
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3.4 Monte Carlo estimators and the radiation field 



3.4.1 Parameterisation of the radiation field 

The treatment of ionization requires tliat tlie radiation field be 
known locally. Since tlie wind may be optically thick, the radia- 
tion iield is computed from the behaviour of the MC packets using 
volume based estimators. In principle, it is possible to describe an 
arbitrary radiation field using MC estimators for the mean intensity 
in narrow frequency bins. However, this requires very large num- 
bers of MC quanta which becomes restrictive for the exploration of 
model parameter spaces. Therefore, we follow previous MC studies 
in adopting a physically motivated parameterisation of the radiation 

field. 

This method has been described bv lLucvl i 19991. 20031 



and used in various other studies (e.g. Mazzali & Lucvl 



20051) 



iLong & Knig"g3l2002 '). In these previous studies, it was assumed 
that the true radiation field was black-body in character and was 
therefore usually parameterised by a radiation temperature and a di- 
lution factor. Here, however, we are concerned with a radiation field 
which is non-thermal in character and so parameterisation based on 
a power-law is adopted: 



no 



2h 



£1 

go Vw^<"^c' 



+ 1 



(15) 



where hvi is the excitation energy of state I above the ion ground 
state. This treatment of excitation avoids the need to specify any 
radiation field parameters beyond those used in the ionization for- 
mula and therefore places no further computational demands on the 
MC simulations. It is very crude but is not of critical importance to 
our study of highly ionized outflows since the excitation state is 
low enough that the ground state populations of H- and He-like 
species are dominant. The most important failing of this treatment 
will be for the metastable triplet states of He-like ions. Thus im- 
provements to the treatment of excitation will be necessary in order 
to extend the code for applicability to less ionized flows and softer 
wavebands. 



3.4.2 Monte Carlo estimators 

The values of W and a are obtained by recording two MC estima- 
tors per grid cell: 



a (r) 



(11) 



where both the parameters W and a are functions of position. 
Given our choice of input radiation field (eguationllOb. if the wind 
is optically thin, we would expect to have a(r) — 1 — F at all 
points. However, for optically thick winds, we expect variation in 
Q. Using this description, the S'i-factors in the ionization equation 
can be expressed as functions of only Te and the local values of W 
and a 



(12) 

This allows the ionization state of the gas to be uniquely de- 
termined from Te, the local mass-density and the radiation field 
parameter pair [W , a] only. 

To specify the excitation state via equation |9] we also need to 
know J I, for ground state transitions. Again, exact MC estimators 
for Jv in all transitions can be constructed but recording and storing 
them for large numbers of transitions in multi-dimensional grids 
becomes expensive in terms of both processor time and memory 
allocation. Ji, can be expressed as 
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(13) 



where cos^^ /i and (p are spherical polar angles, jj,, (j)) is spe- 
cific intensity and /3(^, 0) is the Sobolev escape probability for the 
transition. For the specific case of an homologous flow, (3 becomes 
independent of direction but, in general, it is a function of both /i 
and 0. In the interests of computational expediency, however, we 
neglect this dependency and replace with its angle-averaged 
value, (3. This simplification allows to be expressed as 



dfj. d(^ = /?J, = Wiy°'f3 



(14) 



adopting our parameterisation of the radiation field. This simplifies 
equation|9]to give our final excitation formula 



Ecmf ds 



paths in n 



= ^ EcmfJ^cmf ds 
paths in n 



(16) 



(17) 



where the summation runs over all quanta trajectories which lie in- 
side the grid cell of interest (denoted by 71) and in which the packet 
frequency lies in the regime of interest (Vmi„ < v < fmax). In 
these equations, ds is the length of a Monte Carlo quanta trajec- 
tory length, Ecmf is the comoving frame packet energy and v^mi is 
the comoving frame packet frequency. 

Values for B^^' and E'i'> are recorded during the Monte Carlo 
simulation of all the packets. At the end of the simulation, they are 
used to obtain values of q„ and Wn in each grid cell as follows. 

The ratio £;i^V-En' 

gives the mean frequency. Since we know 
the frequency interval being simulated (fmin < f < I'max), this ratio 
is directly linked to the power-law index via: 



an + 1 



" a„ + 2 
which can be numerically solved for a. 



(18) 



Once a„ is known, Wn can be obtained by normalizing En 



(1). 



Wn = 



E. 



(1) 



(19) 



where Vn is the volume of grid cell n and At is the effective time 
interval represented by the Monte Carlo simulation. 

3.4.3 Iteration 

The formulae given above allow values of Wn and a„ to be ob- 
tained from the behaviour of the MC quanta. However, since the 
values of these parameters also affect the MC simulation (primar- 
ily owing to their importance in determining the ionization state), 
these parameters must be iterated to reach a converged parameteri- 
sation of the radiation field. 
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Table 1. Elements and ions which are included in the radiative transfer cal- 
culations. 



Element 


Ions 


Element 


Ions 


C 


V - VII 


S 


XV - XVII 


N 


VI - VIII 


Ar 


XVII - XIX 





VII - IX 


Ca 


XIX- XXI 


Ne 


IX - XI 


Fe 


XXIII -XXVII 


Mg 


XI - XIII 


Ni 


XXV - XXIX 


Si 


XIII -XV 







Table 2. Inputs parameters for the example model. The upper portion gives 
the adopted physical parameters and the lower portion gives numerical pa- 
rameters. 



In practice, the radiation field parameters are found to con- 
verge very rapidly, often requiring only one iteration step. 



3.5 Atomic data 

Table[T]lists the elements and ionization stages that are included in 
the radiative transfer calculations; for all the light elements, only 
the highest three ionization stages are included while for Fe and Ni 
the five highest stages are treated. 

The solar element abundances of lAsplund et all bOOSi) are 
adopted. 

For all elements other than Fe, atomic models, bound-bound 
oscillator strengths and electron collision rate s were extracted 
from the CHIANTI atomic database, version 5 jPere et alJ[T997l : 
iLandi et al]|2006h . These data include levels with principal quan- 
tum number n ^ 5 for the H-like ions and for states up to the ls5g 
configuration for He-like ions. Ni XXV and XXVI are included only 
for the calculation of the ionization balance; currently no bound- 
bound transitions for these ions are treated. 

For Fe, atomic data were taken from TIPBAs43 which pro- 
vides atomic data from the IRON project for ions of astrophysical 
interest. These data also contain levels for n ^ 5 of Fe XXVI but 
include more highly excited configurations of Fe XXV than CHI- 
ANTI (up to the IslOh configuration). Atomic data for Fe XXIII 
and XXIV were also taken from TIPBASE; however, no states in- 
volving excitation of electrons from the Is shell are included for 
either of these ions. 

Ground configuratio n photoioniza ti on cro ss-sections were de- 



scribed by the fits from Vemer et al. I (Il996l) . except for Ni for 



which fits from I Vemer & YakovlevI ( Il995h were adopted. For ex 
cited states, photoionization cross-sections were computed using a 
hydrogenic approximation for all ions. Since bound-free absorp- 
tion from excited states was not included in the MC simulations, 
these excited state cross-sections are only needed for the computa- 
tion of macro atom transition probabilities (see Section [3.3.4t . The 
fraction of recombinations which go directly to the groun d state , C„ 
was obtained from the hydrogenic calculations of ,Martir> 



4 EXAMPLE CALCULATION 

In this section we present detailed results pertaining to one partic- 
ular instance of our model. 



4.1 Model parameters 

The complete set of adopted input parameters for our example 
model is give in Table (2] 



^ http://vizier.u-strasbg.fr/tipbase/home.html 



Parameter 


Value 


ran Qniirr^p lnminrn;itv f''? 10 l('pV^ T,v 


1 0'^'^ prfTi^ s^-'^ 


Qniirpp nn\x/pr-l iiw nnnfon innpY 1^ 


2.38 


mass of central object, A/i,ji 


4.3 X 10*' Mq 


inner launch radius, r„,i„ 


100 rg = 6.4 X 10^3 cm 


outer launch radius, rmax 


1 ■ S^min 


distance to wind focus, d 




terminal velocity parameter, /i, 


1.0 


velocity scale length, 


''^min 


velocity exponent, /3 


1.0 


launch velocity, no 


0.0 


wind mass-loss rate, M 


0.1 Mq yr-l 


mass-loss exponent, k 


-1.0 


electron temperature, Te 


3 X 10'^ K 


extent of grid 


2 X IQi"' cm 


3D Cartesian RT grid size 


100 X 100 X 100 


2D wind grid 


100 X 100 


Number of MC quanta 


2.4 X 10** 



The black-hole mass, source X-ray luminosity and primary 
power-law index are all chosen to be reasonable for Mrk 766, 
the object for which we make a detailed c omparison in Section [8] 
(A/bh ~ 4.3 X 10*^ Mq |Wangc&Lu"200l']; Lx ~ 10*^ ergs s'^ 
[Poun ds et al]|2003ll ; V ~ 2.38 |Milleretal. 2002]). We assume 
that /u = 1, which is characteristic of radiatively driven flows. 
Furthermore, since we wish to consider outflow features with 
shifts corresponding to v^o ^ 0.1c, we require flow launching radii 
200rg - hence we adopt rmin = lOOrg and rmax = 150rg. 

Physically, i?„ and /3 should depend on the acceleration mech- 
anism and location - however, since such properties are unknown, 
we adopt — rmin and f3 — 1, as appropriate for the simple case 
of an acceleration which occurs on the scale of the system. It is 
assumed that vq <g Voo such that vq — can be adopted in the 
numerical simulations. 

For the example model, we adopt a polar opening angle of 
45 deg (i.e. d — rmin), a moderate mass-loss rate (Al = 0.1 Mq 
yr~^) and electron temperature Te = 3 x 10® K. The effects of 
varying these parameters will be discussed in Section[5] 



4.2 Computed ionization state 

The ionization state is very important since it determines the distri- 
bution of line opacity witiiin the outflow. Fig.[2]shows the computed 
ionization fraction of Fe K-shell ions (i.e. Fe XXV -I- XXVI) in the 
example model. 

There is significant variation in the ionization state in the 
wind. The edge of outflow closest to the rotation (z-) axis is most 
highly ionized since it sees unattenuated X-ray radiation from the 
continuum source. In the outermost parts of the wind the ionization 
state is high and relatively uniform: although far from the X-ray 
source, these regions have low density which disfavours recom- 
bination. Across the wind, there is a significant ionization gradi- 
ent since each layer progressively shields those below it from the 
X-ray source. The very lowest ionization material occurs near the 
xy-plane on the outer edge of the wind - this region has both the 
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Figure 2. The computed ionization fraction for K-sliell Fe (i.e. Fe XXV + XXVI). Only the positive a;2:-plane is shown. Note the logarithmic axes and the 
difference in scale factor on the x- and z-axes. The blue-red colours indicate the ion fraction on a linear scale - red shades represent regions where the wind 
is almost fully ionized (i.e. Fe XXVII is dominant). The black lines indicate lines-of-sight from the X-ray source (the coordinate origin); the lines-of-sight 
drawn are those at the mid-points of the nine angular bins used for the spectra shown in Fig. [5] The solid lines are Hnes-of-sight that never intersect the outflow 
(corresponding to the top three panels in Fig. [5); the dotted Hues intersect the outflow and yield spectra with narrow absorption features (fourth to eighth panels 
in Fig.[3j; the dashed line intersects the outflow and produces a spectrum without narrow absorption features (last panel in Fig. [5). The lines-of-sight which 
intersect the flow are labelled with their integrated Compton optical depth (r) from the continuum source (coordinate origin) to the observer. 



highest densities and the most effective shielding from the X-ray 
source. 

4.3 Computed spectra 

Fig.[3]shows spectra computed for the example model. These were 
obtained by binning the emergent MC quanta by angle relative to 
the polar axis (9). The first eight angular bins shown each cover 
equal solid angle, specifically they encompass A cos 6 = 0.1 while 
the ninth covers the equatorial range < cos 6 < 0.2. 

During the MC simulations, the number of interactions was 
recorded for each quantum. This information has been used to di- 
vide the spectra into "direct" (meaning quanta which underwent no 
interactions) and "scattered" components (note that this includes 
packets that underwent any number of any type of physical event 
- it does not only include true scattering events). These component 
spectra are also plotted in the figure. 

It is apparent that the outflow in the example model affects 
the spectrum as observed from any line-of-sight. The typical ap- 
parent luminosity varies by more than a factor of four - the flux 
level is much higher close to the polar axis than near the equatorial 
plane. This is a consequence of scattering in the flow, primarily by 
free electrons. For small angles relative to the polar axis the source 
radiation escapes directly, unhindered by the wind. At high incli- 
nation angles, the photons enter the wind and most are scattered, 
thereby changing their direction of propagation. Compared to the 
naked source, this enhances the radiation for the polar direction 



and suppresses it close to the equator. Also, as one would expect, 
the scattered radiation field is softer than the direct such that the 
combined spectrum is also slightly softer. 

For viewing angles less than 45 deg the spectrum is unob- 
scured by the outflow and is supplemented by scattered radiation 
which, for this model, can have more than 50 per cent as high a 
flux as the direct component. These lines-of-sight are indicated by 
the solid lines in Fig.|2]and the corresponding spectra are shown in 
the top three panels of Fig. |3] The scattered radiation field shows 
signatures of scattering in spectral lines - a P-Cygni-type line pro- 
files is apparent in the Fe Kq line and there are weaker but similar 
features arising from the Kq lines of lighter elements at softer en- 
ergies. 

At intermediate viewing angles (dotted lines-of-sight in 
Fig. |2j, the line-of-sight intersects the flow. The Compton opti- 
cal depth along these sight lines (r) is of order unity (see labels 
in Fig. [2J meaning that a significant fraction of the continuum ra- 
diation emitted in these directions is scattered by the flow. Dis- 
creet narrow blueshifted absorption lines appear in these spectra 
(fourth to eighth panels of Fig.[3]l. The blue-shift which manifests 
in any particular spectrum depends on the component of outflow 
velocity which lies along the line-of-sight and hence on the inclina- 
tion. In our model, the blueshift increases as 9 decreases. The scat- 
tered component also contains absorption features, albeit weaker 
and broader than those in the direct spectrum. In addition, it in- 
cludes emission which results from line-scattering throughout the 
model and is thus significantly velocity-broadened. 
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Figure 3. Spectra computed for the example model for different lines-of-sight. Each panel shows the spectrum in a different inclination angle bin (the angular 
bin sized are indicated in the figure, measured relative to the polar axis). Note that the ordinate range is different in each panel and does not extend to zero. 
The total spectrum is plotted in black. The spectrum is separated into direct and scattered components and these are overplotted in red and blue, respectively. 
All the spectra are shown normalised to the input primary power-law spectrum. The solid vertical lines indicate the rest energies of Fe Kq for the XXV and 
XXVI ionization stages. 



When viewed from high incMnation angle (dashed line-of- 
sight in Fig. |2] and last panel of Fig. [3), the flow is significantly 
optically thick and the spectrum is totally dominated by scattered 
radiation. These spectra show weak emission lines and weak, rela- 
tively broad absorption lines. 

Thus our example model demonstrates that the simple disk 
wind geometry adopted can naturally produce narrow, blueshifted 
absorption features, as required. However, it also makes clear that 
such features are not necessarily produced and that depending on 
the inclination angle a variety of other features, including relatively 
broad emission or absorption, may be present. 

In the example model, the greatest contrast of the line features 
to the continuum, is only about 40 per cent. There are lines-of-sight 
(e.g. 9 ^15 deg) in which the direct component shows narrow ab- 
sorption features which are close to black at line centre. However, 
in the total spectrum the contrast is always reduced by addition 
of the scattered component in which all features are smeared by 
velocity broadening. We note that our assumption of a stationary, 
smooth flow is likely to maximise this effect since it will generally 
overestimate the relative strength of scattered to direct radiation for 
lines-of-sight in which absorption features occur. 



5 MODEL GRID 

Having discussed the complete spectra from the example model 
(Section|4ll, we now explore the effects of some critical model pa- 



rameters on the computed spectra. To do this we have constructed 
a grid of 45 models and computed spectra for each in 10 viewing- 
angle bins. The parameters which were varied between models, and 
their respective effects on the spectrum, are discussed in the sec- 
tions below. 

In order to elucidate our discussion, we have performed sim- 
ple measurements of the properties of the Fe K absorption features 
in the spectra. First, for each spectrum a continuum spectrum was 
obtained by masking out the strong Ka absorption lines and then 
median filtered to suppress MC noise. The equivalent width (EW) 
of any Fe Kq absorption was then measured by numerical integra- 
tion of the normalised spectrum. The spectra with absorption EWs 
> 40eV were then classified by eye into one of three categories: 

(i) spectra with a single, clear Fe Ka absorption line; 

(ii) those having two, well-separated Fe Ka absorption lines 
(one due to Fe XXV and the other due to Fe XXVI); 

(iii) those with Kq absorption in both Fe XXV and Fe XXVI but 
in which the two lines are strongly blended. 

For the lines in categories (i) and (ii), the blueshift line velocity 
was found by measuring the photon energy of deepest absorption 
and an approximate fuU-width-at-half-maximum (FWHM) was de- 
termined for the absorption feature by fitting a Gaussian profile. For 
the spectra in class (iii), Fe Kq velocities and FWHM were not ex- 
tracted. Although these line properties do not adequately describe 
all properties of the spectra, they provide a convenient quantifica- 
tion of the spectra for the discussions in the sections below. 
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Figure 4. Fe Ka absorption EWs for models with different Af -values 
(all other parameters held fixed). The star, triangle and box symbols in- 
dicate spectra in which there is a single, clear Fe Ka absorption line, 
those in which there is a pair of separated Fe Ka lines and those in which 
the Fe XXVI and XXV components are strongly blended, respectively. The 
colours identify spectra based on the inclination of the observer line-of- 
sight to the polar axis (6). Note that only spectra with Fe Ko absorption 
EWs > 40eV are included. 



5.1 Mass-loss rate, M 

The wind mass-loss rate is the most important model parameter 
which one would wish to constrain observationally. To demonstrate 
its effect, we have computed spectra for models with five M val- 
ues (0.01, 0.03, 0.1, 0.3 and 1.0 Mq yr"^); the Fe Ka absorp- 
tion EWs for these models are shown in Fig. |4] For Mrk 766, 
this range of mass-loss rates explores the region around the mass 
accretion rate estimated from the bolometric luminosity M-^^^ ~ 
Lboi/0.054c^ ~ 0.1 M0 yr~^ a dopting the bolom etric luminosity 
(I/boi = 3 X 10 ergs s~ ) from Hao et al. ('2005^ and the radiative 
effici ency for a Schwarzschild black hole ( Shapiro & Teukolsky 
Il983l) . 

M controls the density in the wind and therefore affects both 
the ionization state and all line or continuum optical depths. At 
the low end of our parameter range, the ionization state is high 
enough that absorption by Fe XXVI is much more important than 
Fe XXV. Increasing M leads to an increase in absorption EWs, as 
one would expect from the increase in opacity. Also for small M, 
strong absorption is only present at large Q since only the very base 
of the flow is dense enough to present significant opacity. 

At intermediate A/, absorption appears across a wider range 
of Q and becomes generally stronger. Increased densities reduce 
the degree of ionization such that absorption by both Fe XXV and 
Fe XXVI Ka becomes more common. Absorption at very high in- 
clination angles becomes less significant since at these angles, the 
spectra are increasingly dominated by scattered radiation which 
washes out any narrow absorption features imprinted on the direct 
radiation. 

For the highest mass-loss rates considered, all the spectra con- 
tain absorption by both Fe XXVI and Fe XXV Ka and scattered 
emission is very important at all inclination angles. This means that 
the EWs are mostly lower in the M = 1 than M = 0.3 yr"^ 
spectra. 



Figure 5. Fe Ka absorption EWs versus line-of-sight inclination for mod- 
els with different d- and M-values (all other parameters held fixed). The 
plotting symbols are as described for Fig. |4] Note that only specfi'a with 
Fe Ka absorption EWs > 40eV are included. 
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Figure 6. The blueshift of the Fe XXVI Ka line versus line-of-sight in- 
clination for models with different d- and Af -values. The plotting symbols 
and colors are as described for Fig [5] Note that only spectra with Fe Ka 
absorption EWs > 40eV are included. 

5.2 Wind geometry, d 

Within our simple wind prescription, d controls the opening angles 
of the outflow and therefore determines both the angular range of 
lines-of-sight which pass through the wind and the relative orien- 
tations of the outflow and rotational velocity fields. To explore the 
effect of geometry on the spectrum, we have extended our grid of 
models to include d-values of 0.33 rmin and 3 r-min, thereby produc- 
ing geometries which are, respectively, more equatorial and more 
polar than the example model. For these d-values, we have com- 
puted models with the same set of A/- values used in Section [STI 
keeping all other parameters fixed to those of the example model. 

Fe Ka absorption EWs are plotted for a sample of these mod- 
els in Fig. [5] This shows the role of the wind opening angles in 
determining which lines-of-sight intersect the outflow in regions 
where optical depths are high enough to cause line absorption. 

The wind opening angles and line-of-sight orientation also 
affect the observed Fe Ka blueshift. Fig. [6] shows the Fe XXVI 
blueshift velocity for the same models used in Fig. |5] excluding 
spectra with blended Fe XXV and XXVI Ka absorption. As men- 
tioned in Section |4l the variation in blueshift arises from line-of- 
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Figure 7. Fe Ka absorption EWs versus line-of-sight inclination for mod- 
els with different r,^i„- and M -value. The plotting symbols are as described 
for Fig.|4] Note that only spectra with Fe Ko absorption EWs > 40eV are 
included. 



Figure 8. Fe Ka absorption EWs for models with different Tg-values (all 
other parameters held fixed). The plotting symbols are as described for 
Fig|4] Note that only spectra with Fe Ka absorption EWs > 40eV are 
included. 



sight projection of the outflow velocity and the effect is large: there 
can be more than a factor of two spread in the apparent blueshift. 
Since the maximum outflow velocity is 



2GA/bh 



= 0.14c = 4.2 X lO^km s"^ 



(20) 



for all the models discussed here, the variation in blueshift is 
entirely a consequence of orientation. Thus, while the apparent 
blueshift can be as large as the terminal outflow velocity, it is often 
significantly smaller. 



5.3 Inner launch radius, rmi„ 

The inner launch radius of the wind, r^m, determines how far the 
wind material is from the X-ray source and affects the wind den- 
sity since the surface area of annuli in the disk plane increases with 
r. To investigate its effect, we have computed spectra for models 
with r„,in larger than in the standard model {r^m = 200rg and 
500rg) adopting the same range of i\/-values used in Section [STI 
For all these models, we also varied other parameters as neces- 
sary to preserve the relationships used for the standard model (i.e. 
r-raax = l-5rmin; d = rn,i„; Rv = Tn^m, fv = !)• Fig.|7]shows Fc Ka 
absorption EWs computed as a function of viewing angle for six of 
these models. 

Changing rm[„ affects the typical range of viewing angle for 
which absorption features are strong, favouring smaller cos 9 for 
larger rmin- For fixed M, increasing rmin tends to make the absorp- 
tion lines weaker, a consequence of the reduced density in the flow. 



5.4 Electron temperature, Te 

Under our assumptions, the only important role played by the elec- 
tron temperature is in determining the ionization fractions. To il- 
lustrate its effect, we have used models which span the same M- 
grid as discussed in the section above but which adopt Te both 
lower (lO'' K) and higher (10^ K) than used in the example model. 
EWs for a subset of these models, namely all those having M = 
0.1 Mq yi~^ are shown in Fig.[8] 

The dependence on Te is these models is very simple. As Te 
increases, the wind becomes more ionized, reducing the H- and 
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Figure 9. Fe Ko absorption EWs for models with different F-values (all 
other parameters held fixed). The plotting symbols are as described for 
Fig|4] Note that only specfi'a with Fe Ka absorption EWs > 40eV are 
included. 



He-like ion populations in favour of the fully ionized state. This re- 
duces the line opacities and causes the absorption EWs to decrease. 
Absorption in both Fe XXV and XXVI is most common at low Te 
since the Fe XXV population is depleted more rapidly than that of 
Fe XXVI as the wind becomes more ionized. 



5.5 Primary power-law index, F 

The source power-law index determines the fraction of the X-ray 
flux which is hard enough to contribute to the ionization of Fe XXV 
and XXVI and therefore has an effect on the Fe Ka absorption EWs 
which is qualitatively similar to that of Te . In Fig.|9]we show Fe Ka 
EWs for models with M = 0.1 Mq yr"^ adopting F = 1.88, 2.38 
and 2.88, with other parameters as in the example model. As ex- 
pected, higher Fs lead to larger EWs since there are fewer Fe K- 
shell ionizing photons and hence the ionization state is lower. 

Note that the total 2-10 keV luminosity (Lx) is constant here 
- although we have not explicitly explored the effect of changing 
Lx, it is clear that it will have a similarly simple effect on the 
ionization state. 
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Figure 10. Fe Ka absorption EWs versus the Fe XXVI blueshift. All spec- 
tra computed from our grid of models (450 specti'a in total) with unblended 
Fe Ka absorption with EW > 40 eV are included in the plot. The star 
symbols represent model spectra in which a single Fe Ka absorption line is 
dominant. Triangles indicate spectra with two, unblended Fe Ka absorption 
lines. The solid lines coimect points for the same model but viewed from 
different inclination a ngular bins. The re d cross indicates the Fe Ka absorp- 
tion line reported by Pounds et al. (2003) for PG1211+143 while the red 
bars indicat e the range of EWs obtained with the t ime sliced observa tions 
of Mrk 766 ^Turner et alj2007l) and MCG-5-23-16 Ifiraito et al.l2007l) . 



6 OBSERVABLE PROPERTIES 

In Section[5]we constructed a grid of models to explore the effects 
of various wind properties on the Fe Ka absorption line properties. 
However, since the relevant wind parameters are not known a pri- 
ori it is instructive to investigate the properties of the models in the 
space defined by observable quantities. Given the typical data qual- 
ity and limited spectral resolution of current X-ray telescopes, the 
primary observables pertaining to Ka absorption lines are the EW 
and the velocity shift of the lines. 

Of the 450 spectra which were generated with our grid of mod- 
els (Section ID, 101 were classified as having Fe Ka absorption 
EWs of > 40 eV and, of these, 79 were classified as having ei- 
ther a single or an unblended pair of absorption features such that 
a meaningful velocity blueshift could be extracted from the spec- 
trum. Fig.[TO]shows the distribution of Fe Ka absorption EWs ver- 
sus the Fe XXVI blueshift for these 79 spectra. 

There are a wide variety of Ka absorption features in our mod- 
els. EWs of more than 100 eV are rare but do occur; the maximum 
EW produced by any of the models is ~ 135 eV. As one would 
expect, the EW tends to be larger when two lines are present al- 
though there are cases where a single line achieves an EW of nearly 
100 eV. As already highlighted in Section [J!2l a very wide range 
of blueshift velocity is predicted. Typically, the maximum EWs is 
larger when the blueshift is larger; but there is significant variation 
in EW at all blueshifts. In general, there are many sets of wind pa- 
rameters that can produce similar EWs and blueshifts such that Ka 
properties cannot uniquely discriminate between our models. Other 
potential outflow signatures are discussed in Section|7] 

As a simple comparison, points representing reported de- 
tections of narrow blueshifted Ka absorption in the spectra of 
three real AGN (PG121 1-1-143, Mrk 766 and MCG-5-23-16; see 
below) are overplotted in Fig. [TO] Observations clearly show 
that outflow absorption f eatures are signifi c antly time variable 
(e.g. iRisaliti et al.l ]2005l ; iMiller et al.l |2007| ; iTutrier et alj |2007| : 
iBraito et ak 2007 ) and thus one must be cautious when compar- 



ing them directly to the results of our steady-state models - the de- 
velopment of multi-dimensional models involving clumpy or struc- 
tured flows is necessary for a detailed comparison. 

Pounds et al. (200 3i) first reported blueshifted Ka absorption 
from H- or He-like Fe in the spectrum of PG121 1-1-143 (see also 
[Pounds et al]|2005l; [Pounds & Pag3l2006l ; [Pounds & Reevesll2007l 
for further discussion of this interpretation of the PG121 1-1-143 
spectrum); they rep ort an absorption EW of 105 ± 35 eV 
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2005h and a blueshift of either 0.08 - 0. 10 c 
l2005l) or 0.13 - 0.15 c jPounds & Pag3l2006h . de- 



pending on whether one assumes the observed features is owing 
to Fe XXVI or XXV. The point plotted in Fig. [lO] for this object 
assumes the Fe XXVI interpretation a nd the measurement uncer- 
tainties quoted by [Pounds et alj ( [2005[) . 

Two distinct absorption feat ures (6.9 and 7.2 k e V) were re- 
porte d in spectra of Mrk 766 jTumer et all [2007[ ; [Miller et all 
[2007h . at 6.9 and 7.2 keV. Although the nature and relationship 
between these features is unclear, one possible interpretation is that 
they correspond to blueshifted absorption by Fe XXV and Fe XXVI, 
respectively - in t his case, the blueshift velocity of both lines is 
~ 13, 000 km s"^ jMiller et al.l2007h . [Tumer et al.[ ( l2007h showed 
that the strength of these features varies significantly in time. The 
largest fluxes of the absorption features occurred in their time slice 
9 (see their fig. 2) and had an EW ~ 80 eV each. In some time 
slices, the line fluxes are consistent with zero. In Fig.[TOl the range 
of Mrk 766 absorption feature properties is marked adopting the 
13,000 km s^^-outflow interpretation (i.e. identifying the 7.2 keV 
feature with Fe XXVI) and the EW range of - 80 eV as motivated 
by t he time-sliced data. 

[Braito et"ai] l [2007[) detected a variable absorption feature at 
7.7 keV in the spectrum of MCG-5-23-16. This feature is also indi- 
cated in Fig. [To] adopting the upper bound of on its EW range based 
on their analysis of the one of five 20 ks time intervals in which the 
line was most clearly detected (see their table 4). The velocity shift 
was computed assuming that Fe XXVI is dominant. 

In terms of shift and range of EW, the absorption features in 
PG1211-F143, Mrk 766 and MCG-5-23-16 fie close to the space 
explored by our grid of models (Fig.llOll. Thus wind models such 
as those we have investigated are promising starting points for ex- 
plaining the Ka absorption in these observations. 

Blues hifted Kg absorptio n features have also been reported in 
NGC 1365 jRisaliti et al.l2005h . In this case, the shifts are consider- 
ably smaller, ~5000 km s~^ at most. Although our grid of models 
does include some spectra with features at these low velocities (see 
Fig. (To), the EWs in the computed spectra are below those observed 
by a significant f actor (the combined Fe Ka EWs reported by 
[Risaliti et alj2005l are in the range 200 - 300 eV, lying off the scale 
of Fig. not . Even lower-blueshi ft highly-ionized absorption has 
been i -eported in MCG- 6-30-15 jYoung et a l 2005; Miniu tti et al.l 
[20071 ; [Miller etalj[2008[) and NGC 3516 frurner et al..,200i) . The 
simplest way in which such strong, relatively low-velocity features 
might be produced with our model would be by considering lower 
terminal flow velocities but it may also be possible to obtain deeper 
absorption by considering flows which are significantly clumped. 
Detailed study of these possibilities lies beyond the scope of this 
paper, but we do note that even the NGC 1365 features are within a 
factor of 5 in EW of those obtained with our grid of models. 

In principle, the model spectra indicate that the absorption line 
shapes are complex and therefore provide information on the out- 
flow properties. However, given the limited sensitivity and spectral 
resolution of current X-ray telescopes, determining line shapes is 
virtually impossible - at best, only a rough limit on the characteris- 
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Figure 11. Fe Ka FWHM versus the Fe XXVI blueshift. The plot symbols 
are the same as in Fig. llOl 



the spectra computed from our grid of models, particularly those 
with high A/-values. The strongest of these Ka features are those 
of S, Si and O - the absorption EWs for these lines are up to 10, 10 
and 3 eV respectively. Their strengths correlate closely with each 
other and more loosely with the Fe Ka feature. 

At energies harder than Fe Ka, absorption by H-ZHe-like Ni 
Ka and by Fe K/3/K7 occur. Individually, these features are all 
weaker than Fe Ka but they are often blended together producing 
fairly significant absorption which can extend up to around 10 keV. 

These additional absorption lines can provide both confirma- 
tion of outflows detected via Fe Ka and supplementary constraints 
on the flow properties. However, their reliability is limited since (i) 
in our models, they are always weaker and thus harder to detect 
than Fe Ka and (ii) with the exception of the hard energy Ni or 
Fe K/3 lines, their identifications are less clear since there are more 
potential atomic transitions at soft energy. Thus, as diagnostics for 
highly ionized flows, these lines are always of subordinate value to 
FeKa. 



tic linewidth is available. For comparison with such limits. Fig. II II 
shows the FWHM determined from our spectra (by Gaussian fit- 
ting) versus the Fe XXVI blueshift for the same models as plotted 
in Fig. [To] 

The FWHM varies between about 3.5 x 10^ and 1.8 x 
lO* km s~^, with values around 5-7 xlO^ km s^^ being most 
common. For all our model spectra in which the blueshift is greater 
than 10^ km s^^, the FWHM is smaller than the blueshift mean- 
ing that the absorption does not extend down to the rest energy 
of the line. As one would expect, the largest linewidth (~ 1.8 x 
10* km s"^) occurs in model where d has the smallest value con- 
sidered (d = 0.33) and for a viewing direction which has a large 
pathlength through the flow. 

To our knowledge, there are currently no firm measurements 
of linewidths for narrow bluesh ifted Fe Ka a bsorption features in 
AGN spectra. For PG121 l-l-143. |Pounds et afl {2003) placed a limit 
of < 1.2 x 10* km s~^ on the Fe Ka linewidth, which is consistent 
with th e majority of the m easurements from our spectra shown in 
Fig. fmlBraito et alj ( l2007l) found a best-fit line width of cr = 200± 
100 eV (FWHM = 1.8 x 10* km s"*) for the absorption feature in 
MCG-5-23-16. However, the fit in which they allowed the linewidth 
to vary was not a significant improvement over that in which they 
piimed the width at a = 100 eV. Thus, their line width constraints 
are also consistent with the majority of the spectra described by 
Fig.E] 



7 OTHER SIGNATURES OF A FAST OUTFLOW 

As discussed in Section [T] blueshifted Fe Ka absorption is the 
clearest signature of highly ionized outflowing gas and has there- 
fore been the central topic of this paper. However, as mentioned in 
Section|4l narrow blueshifted Fe Ka absorption is neither a neces- 
sary consequence nor the only signature of a highly ionized flow. 
Thus if such outflows are common in the AGN population it is im- 
portant to consider other effects they may have on X-ray spectra. 
In this section, we briefly describe some of these effects and their 
possible implications with reference to our model spectra. 

7.1 Blueshifted absorption lines other than Fe Ka 

Although the Fe a feature is typically the most prominent absorp- 
tion line, Ka absorption by lighter elements is present in many of 



7.2 Emission lines and P-Cygni profiles 

As pointed out in the reference spectra shown in Section |4l the 
blueshifted absorption lines are often accompanied by broad emis- 
sion which is centred close to the line rest energy. Since these emis- 
sion features are formed as a result of line scattering and recom- 
bination emission in the outflow, their strength strongly depends 
on the density of the wind; thus although they are fairly weak in 
the spectra from the example model, they become more signifi- 
cant in our models with higher mass-loss rates. Reliable predic- 
tions for these phenomena in non-spherical outflows requires multi- 
dimensional radiative transfer and thus our methods are particularly 
well-suited to describing them. 

To illustrate the types of line profile, Fig.ll2lshows a montage 
of spectra from models with relatively high mass-loss rates. There 
is a wide range of Ka profile types in these models and their emis- 
sion EWs, measured relative to a power-law continuum fit, can be 
in excess of 200 eV; this is compara ble to the typical bro ad Ka 
EWs measured in a sample of AGN bv lNandra et al] ( |2007|) . 

At low inclination angles, our model Fe Ka profiles have a 
near-P-Cygni character: they shows some blueshifted absorption 
and broad, redshifted emission, the strength of which grows with 
M. Such profiles a re qualitatively sim ilar to those observed in some 
Seyfert 1 galaxies dOone et al.]|2007h . 

At intermediate angles, when the line-of-sight is close to look- 
ing down the wind cone, the narrow absorption lines are strongest 
but are accompanied by moderately strong, emission lines with 
fairly extended red wings. We note that the narrow absorption lines 
become less prominent at very high A/, a consequence of the in- 
creased contribution of scattered radiation in the spectra. 

For the highest inclinations and Af-values, the Fe Ka line is 
almost pure emission, peaking at the rest energy of the line and 
having a long red tail. The redward extension of the line emission is 
not a consequence of gravitational redshift but arises from electron 
scattering in the outflow. 

Although somewhat more diverse, our red-skewed emis- 
sion line profile shapes are qualitatively similar to those ob- 
tained via the s ame physical mechanism in spherical outflow 
models (i Laurent & Titarchuk 2007; see also iLaming & TitarchukI 
[2004) whil e they are significantly broader than those computed 
by .Rozanska & Madeil ( [2OO81) for Compton scattering in irradi- 
ated accretion disks. The la rgest Fe Ka emission EWs found by 
iLaming & TitarchukI ( |2004|) ( ^ 4keV) substantially exceed those 
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Figure 12. Sample spectra computed for different viewing angle bins (left to right) for models with differing mass-loss rates (top to bottom). Except for M, 
all model parameters are fixed at those of the example model. The plotted spectra are all normalised to the input primary power-law spectrum. 



found in any of our models but we note that the large EWs in their 
models mostly arise from significantly lower ionization stages than 
are present in our models. 

These emission line properties show that highly ionized out- 
flow may affect the Fe K region beyond imprinting narrow, 
blueshifted absorption lines. This may have consequences for study 
of the Fe K fluorescent emission which originates in AG N accre- 
tion disks (e.g. | Fabian et al. 1989: Nandra et al. 1997; Fa bian et al.l 
l200d : lMineill2007l : iNandra et alj|2007l) . since it may contaminate 
the disk emission and /or lead to an apparently multi-component 
emission line (see e.g. lO'Neill et al.l 2007). In this context, the red 
wings predicted for the emission lines may be of particular rel- 
evance in view of the potential for confusion with the effects of 
gravitational redshift ~ however, more sophisticated 3D models go- 
ing beyond the smooth-flow assumption of our parameterised wind 
models must be examined before such a possibility can be consid- 
ered in much greater detail. 



7.3 Spectral curvature 

Absorption by outflowing material has been discussed as a 

possible explanation for the so -called soft-excess in X-ray 

spectra {o ierlihski & Done 200l l2006l : iMiddleton et"ai] 120071 : 
ISchurch & Donell2007l . I2008IV In this picture, the decrease in flux 
above ~ 1 keV is attributed to absorption by light or intermediate 
mass elements. 

Such absorption does occur in our models, particularly for 
high M-values (see Fig.ll2t however the scale of the effect is too 
small in the models presented here: the typical observed soft excess 



requires a drop in flux of nearly a factor of two between about 1 
and 2 keV jMiddleton et al.ll2007h . Furthermore. ISchurch & Done! 
(2008) have argued that very large, relativistic velocities are re- 
quired for the absorption model to work since the observed soft 
excess spectra are very smooth. Our current models support their 
conclusions since all cases in which significant spectral curvature 
arises are accompanied by discreet features. 

Further consideration must await extension of our calculations 
to low-ionization states since more significant absorption at rela- 
tively soft energies can then be expected. 



8 DETAILED COMPARISON WITH MRK 766 

One of the AGN whose X-ray spectrum pro vided some of the 
initial motivation for this work is Mrk766 jMiller et al.l l2007l : 
iTumer et alj2007h . As well as t he possible outflow absorption lines 
already mentioned in section|6l iMiller et"ail ( l2007t) suggested that 
the extremely hard 2-10 keV spectrum in the low flux state, cou- 
pled with the strong Fe edge and weak Fe emission line, might be 
a signature of either ionised absorption or ionised reflection, pos- 
sibly from an obscuring wind. Such a picture provides an alterna- 
tive explanation to the hypothesis of "relativistic blurring" that has 
been invoked to explain the apparent "red wing" at energies below 
7 keV in the mean X-ray spectra of AGN, including Mrk766 (e.g. 
iNandra et alj2007l) . 

We now make a direct comparison between the wind model 
in this paper and the mean spectrum of Mrk 766 obtained by sum- 
ming the XMM-Newton EPIC PN data from the 2000, 2001 and 2006 
observations. Observation IDs were 0096020101, 0096020101 
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current model can explain the full spectral variability seen in this 
AGN's various flux states: the hard spectrum in the low state is most 
likely caused by m oderately ionised material, ^ ~ 100 erg cm s~^ 
jMiller et alj2007h . which is not yet treated in the wind model. We 
postpone to future work fits to the X-ray spectra of this and other 
AGN with enhanced, lower-ionization, wind models. 



/keV 



Figure 13. The best-fitting wind model fit to the mean 2-10 keV spectrum 
of Mrk766, plotted in units of Ef(E), showing model and "unfolded" data. 



and 0304030[1-7]01 and data were processed as described by 
jMiiler et ala,2007.) . As in that work, energy bins of width equal to 
the energy-dependent HWHM were adopted. The model was fitted 
to the data by first creating a multiplicative table ("mtable ") cover- 
ing a range of parameters that could be used with XSPEC jArnaud 
[199^ . The models used for this table were those described in Sec- 
tions 5.1, 5.2 and 5.3 making up a regular, four-dimensional grid 
covering three of our most important wind parameters (rmin, d and 
M) and the viewing-angle (6). The other wind parameters were not 
varied in the fit. The table was used to multiply an assumed in- 
trinsic power-law whose photon index was a free parameter in the 
fit: in principle the fitting is valid only for the choice of power- 
law index used to create the model, but in practice we find that 
there is little variation between models whose power-law indices 
differ, and that the variation that is seen is to some extent degen- 
erate with the effects caused by varying other wind parameters 
(see Section[5](. We expect also some amount of spectral curvature 
caused by the ionised absorption previously identified in Mrk 766, 
so w e add a low column of ionised absorption calculated by XSTAR 
jKallm an et al. 2004.) and we model the narrow 6.4 keV Fe emis- 
sion line as arisi ng in mildly ionised re flection calculated by the RE- 
FLIONX model ( IRoss & Fabianll200^ with no relativistic blurring. 
Owing to computing time limitations the model grid was necessar- 
ily rather coarse in parameter space, so to check the fit of the best- 
fitting model, interpolated between models, we then recalculated 
the model with precisely the parameter values indicated by fitting. 
The wind best-fit parameter were rmin = 385 rg, d = l.Brmin, 
M = O.4M0yr"^ and 008 6* = 0.77, while the best-fit pho- 
ton index was F = 2.03 and absorber parameters were Nh = 
5.3 X 10^^ cm~^ and ionization parameter ^ — 40ergcms~^. 
These parameters yielded a goodness-of-fit ~ 145 for 108 de- 
grees of freedom (no component of systematic error is assumed). 
The fit to the spectrum in the 2-10 keV range is shown in Fig.ll3l 
where we plot the model spectrum at the XMM-Newton resolution 
together with "unfolded" data points. 

The fit of the wind model to the data does indicate that this 
may well provide a good explanation of the shape of X-ray spectra 
in AGN such as Mrk 766. With the exception of the narrow Fe emis- 
sion line, most of the spectral shape in the Fe regime is provided by 
the wind model, with the reflection component contributing only 
3 per cent of the total 2-10 keV energy flux and only 4 per cent of 
the flux density in the continuum around Fe. It is unlikely that the 



9 CONCLUSIONS AND FUTURE WORK 

We have used a simply-parameterised model to construct a variety 
of possible AGN outflow geometries and computed spectra for vari- 
ous lines-of-sight using a multi-dimensional Monte Carlo radiation 
transport code. 

Our results show that narrow blueshifted Fe Ka absorption 
lines are a natural consequence of highly ionized flows and that 
their formation depends strongly on the flow inclination angles rel- 
ative to the observer's line-of-sight. Typically, our models suggest 
that moderately strong features Ka absorption features ( ^ 40 eV) 
will manifest for a significant minority of possible viewing angles. 

Although blueshifted absorption lines are the clearest signa- 
ture of an outflow, our calculations have illustrated a number of 
other spectroscopic signatures which may be relevant to the inter- 
pretation of X-ray data in the context of AGN outflows. In particu- 
lar, electron scattering, line scattering and recombination emission 
can all affect the Fe K region. For suitable wind conditions and 
lines-of-sight, these processes can lead to significant net emission 
in Ka and can make a strong red electron-scattering line wing. We 
have illustrated that these emission line profiles look qualitatively 
similar to observed spectral features (e.g. in Mrk 766) suggesting 
that outflows may have a much more general role in the formation 
of AGN X-ray spectra than simply imprinting absorption lines. At 
the very least, emission/scattering features such as those found in 
our models must be considered as a potential source of contamina- 
tion for studies of other sources of Fe K emission. 

Considerable further work is required to more fully under- 
stand the role of outflows in the spectral formation of AGN. The 
primary limitation of this study is the restriction to reliable mod- 
elling of only K-shell ions. K-shell transitions occurring in lower 
ions, as may be present in less ionized wind regions, will have a 
qualitative effect on the Fe K region - e.g. while the Ka line of 
H- and He-like ions predominantly scatter incoming photons, the 
Li- and Be-like ions can destroy them with significant probability. 
Fortunately, the restriction to K-shell ions is not a fundamental lim- 
itation of the methods used but was made primarily in the interests 
of computational expediency for this study of Fe K shell absorp- 
tion features. The next step will be to extend the physics of the 
code to incorporate L-shell ions so that lower ionization states may 
be probed. In particular, this will allow the investigation of outflow 
geometries with larger r„,ax and flows with a clumpy structure. Fur- 
thermore, in this study we have only considered the simplest X-ray 
source location and have neglected any interplay between differ- 
ent components in the AGN environment. In particular, we have 
not considered how emission, absorption and reflection by an ac- 
cretion disk will affect the spectral signatures of outflow, and vice 
versa. Ultimately, it will be desirable to construct physically mo- 
tivated models for complete disk/outflow systems and use these to 
interpret the increasing wealth of X-ray observations. 
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